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1. Introduction 

Interest in precise lattice calculations of K — > %% decays stems from the possibility of gaining 
insight into the origin of the AI = 1/2 rule, and of putting constraints on CKM matrix elements. In 
particular, e'/ e and the CP violating phase can be computed when lattice calculations are compared 
with experimental results [jl], |[] . 

Several quenched K — > 7T7T calculations have been performed in the past [^|, ||, [| ||, §]. The 
calculation presented here is a quenched calculation using domain wall fermions (DWF) on a 
24 3 x 64, = 16 lattice with strong coupling (a -1 = 1.3 GeV). This represents a larger spa- 
tial volume, (3.6 fm) 3 , than in previous works, allowing for a lighter, near physical pion mass of 
m % = 227.6(6) MeV. As in Refs. [^, f|, ^ |6|] we are directly calculating the two pion decay ampli- 
tudes and not relying on chiral perturbation theory to determine these amplitudes from the simpler 
K — > % and K — ► |0) amplitudes. This use of chiral perturbation theory, upon which Refs. [[I], 
were based, now appears to introduce large, uncontrolled systematic errors [pi]. This quenched 
study is a pilot project that will be followed by a full QCD calculation using the RBC/UKQCD 
32 3 x 64, L s = 32 lattices with 2+1 flavors of DWF and strong coupling (a -1 = 1.4 GeV) that are 
currently being generated, with light pion masses (m n < 250 MeV). 

2. Four Quark Operators and the Effective Hamiltonian 

The weak interactions and the effects of the heavier quarks can be included in the lattice QCD 
simulation by evaluating matrix elements of an effective Hamiltonian [|8[ [J. In particular we use 
the conventions of equation 3 in |0|. In this formalism we are interested in calculating matrix 
elements of four quark operators between K and nn states. Here we calculate matrix elements 
with AI = 3/2 operators, where AI is the change in isospin. The AI = 3/2 operators are further 
classified by how they transform under SU(3) L x SU(3) R into "(27,1)", "(8,8)" and "(8,8) mixed" 
parts. [|,|]. 

3. Twisted Boundary Conditions 

The kinematics of the physical K — > %% decay are such that the pions have significant mo- 
mentum, ~ 200 MeV each, so that it is required to simulate a two pion excited state on the lattice. 
Extracting an excited state requires multi-exponential fitting making it difficult to obtain a signal. 
One remedy is to use twisted boundary conditions [|[ f|, ||J [l0[ ], in which the boundary conditions 
on a quark field, instead of simply being periodic, cause it to change by a phase e"^ when going 
through the boundary. We say that this spatial direction is "twisted" by an amount (p. 

Twisting a quark by (j> will cause it to have momentum 

Pn = £ ( 3 -^ 

where L is the spatial extent of the lattice. In particular, we use twists of amount % (corresponding 
to antiperiodic boundary conditions in the twisted direction) to obtain a two pion system with 
total momentum. For example, if only the x direction is twisted then we can have a two pion ground 
state in which one pion has momentum p x = ? and the other pion has momentum p x = — 5 . We 
can obtain two pion ground states with pion 3-momentum of magnitude 

n V2n V3n 

p * = i — — (3 - 2) 
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by twisting one, two, and three spatial directions respectively. Since we actually perform the cal- 
culation for the decay K + — ► n + 7i + , and relate it to the physical decay via the Wigner-Eckhart 
theorem [^J|, we can twist the down quark only, thus giving the two pions momentum and giving no 
momentum to the kaon. By rotational symmetry we can average over results for pion momenta of 
the same magnitude but different directions, i.e. we can average over results with the same number 
of twists. 

4. Details of the Calculation 

The calculation was carried out on 68 configurations of quenched 24 3 x 64 lattices using 
the DBW2 action, and domain wall fermions with L s = 16. We used the QCDOC computers 
at Brookhaven and Columbia University for both configuration generation and propagator inver- 
sion. The inverse lattice spacing is a -1 = 1.3 GeV, the physical volume is (3.6 fm) 3 , and we set 
the light and strange quark masses to m/ = 0.0055 (chosen so that m K L « 4) and m s = 0.08 in 
lattice units respectively. This yielded a pion mass of m K = 227.6(6) MeV and a kaon mass of 
rtiK = 564(2) MeV in physical units. The pion momentum corresponding to one twist for these 
lattices is p n = 7i/L = 170 MeV so that the decay is expected to be nearly energy conserving 
(E„ n -m K & 2 y /m£+pj'- m K = 4 MeV). 

We added and subtracted quark propagators with periodic and antiperiodic boundary condi- 
tions in the time direction from each other in order to double the effective time length. Periodic 
plus antiperiodic (P+A) propagators thus provide a source at t = and periodic minus antiperiodic 
(P-A) propagators provide an effective source at t = 64. These provide the left and right walls for 
the kaon and two pions and the time t of the four quark operator is varied. The kaon mass, two 
pion energy, and K ^ KK matrix element were extracted from the asymptotic behavior of the cor- 



responding correlation functions as in [11], the only difference being that the K — » %% correlator 



was fit directly rather than fitting the quotient of correlators defined in equation 3.5 of Jl lfj - 

For the kaon and the two pions with zero momentum we used propagators with Coulomb 
gauge fixed wall sources. For the two pions with non-zero momentum, we used the same type 
of propagators for the u quark, but used propagators with twisted (antiperiodic spatial) boundary 
conditions for the d quark with Coulomb gauge fixed momentum wall sources of the "cosine" type 

Sp,co S ine(x) = COS(/7^x) COS(p y y) COs(p z z) (4. 1) 

Had we used sources of the "pure momentum" type 

s p (x) = e* x . (4.2) 

it would have been necessary to do an inversion with source momentum +p for one d quark and 
an inversion with source momentum — p for the other d quark in order to couple to a two pion state 
with zero total momentum. Using the same cosine source propagator for each d quark eliminates 
the need for this extra inversion, but this leads to cross terms that couple to two pion states with 
non-zero total momentum. For example, with a single twist in the x direction, we hope to couple 
to a two pion state with individual pion momenta pi = — fx and P2 = fx, and thus the product of 



the sources of the two d quarks according to Equation 4. 1 should be 

V, cosine (xi>p 2 ,cosine(x2) = COS (j^Xl) COS (^-X 2 

1 



_ ( e 'i A Vi* 2 +e'z Xl e- l * X2 +*>- I r J V^ 2 + e^i* 1 e^* 1 ) (4.3) 
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Table 1: Lattice matrix elements for a time separation of 64 between the kaons and the two pions (first 
calculation). Note that these values differ from those presented in the talk where a normalization factor was 
omitted. 



p 


^(27,1) 


^(8,8) 


•^(8,8) mixed 




k/L 


0.00408(40) 


0.0950(77) 
0.161(85) 


0.294(34) 
0.55(30) 



But in Equation 4.3 the first and last term couple to two pion states with total momentum 2 j- and 
—2 j- in the x direction respectively. To circumvent this difficulty, we calculate two pion correlators 
with cosine sources but pure momentum sinks. The unwanted terms from the source vanish in the 
limit of an infinite number of configurations since the two pion final state is constrained to have 
zero total momentum by the pure momentum sink, but since these are sinks they require no extra 
inversions. In the K — > 7171 correlator, the momentum kaon in the initial state has a similar effect 
on the cosine sources of the two pions in the final state. 

5. Results of the First Calculation 

The calculation was carried out as described, and the results for the matrix elements are shown 
m Table[I[ The 

errors on the (8,8) and (8,8) mixed matrix elements with momentum 7t/L (where 
L = 24) are of order 50%, and the data was too noisy to extract a signal for the (27,1) matrix element 
with momentum n/L and for all the matrix elements with momentum \fl%jL and s/3n/L. This is 
due to the fact that a time separation of 64 is too large, especially for the higher momenta where 
the higher energy of the two pion state causes its signal to decay into noise well before it can reach 
the four quark operator. 

6. Results of the Second Calculation 

The calculation was redone with smaller time separations between the kaon and the two pions. 1 
Specifically, propagators with sources at tg =20, 24, 28, 32, 36, 40, and 44 were generated in order 
to calculate K — > 7T7T correlators with kaon sources at these times, while the two pion sources 
remained at either the left or the right walls. Thus we could achieve time separations between the 
kaon and the two pions of 20, 24, 28, and 32 in two different ways in order to double the statistics 
(for example a time separation of 20 results when the two pions are at the left wall, i.e. t K = 0, 
and when tx = 20, and also results when the two pions are at the right wall, i.e. t„ = 64, and 
when tg = 44). These separations were chosen so that the kaon and two pion sources were close 
enough to the four quark operator (with the operator between the sources) so that the signals from 
the sources do not decay into noise before reaching the operator, but were sufficiently distant from 
each other that a plateau region can be found where the effects of possible excited states have died 
away. 

Effective mass plots for the kaon correlator, the two pion correlator with momentum 7i/L, and 
the K — > Tin correlator with the (27,1) operator, momentum 7i/L and time separation 24 between 
the two sources, are shown in Figures [I] and ^[ The two pion energies that were extracted from 
these effective mass plots are given in Table ||. Table || shows the lattice matrix elements obtained 



These results were not presented in the talk. 
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for different time separations of the two sources, as well as an error weighted average of the results 
from the four different time separations. The error of the error weighted average is smaller than the 
error of the matrix element for any single time separation in all cases. For momentum and n/L 
the error of the error weighted average is of order 3-4% (for all of the operators), for momentum 
y/2n/L it is of order 7%, and for momentum yhizjL it is of order 15%. 




Figure 1: Effective mass plots for the kaon correlator (left) and the two pion correlator for pions with 
momentum n/L (right). The effective mass of a correlator C(t) is defined as m e f[(t) = — In (C(t)/C{t — 1)). 




Figure 2: Effective mass plot for the K — > nn correlator with the (27,1) four quark operator, pions with 
momentum n/L, the two pion source located at t = and the kaon source located at t = 24. 

7. Conclusion 

Lattice matrix elements were computed for AI = 3/2, K — > 7171 decays on quenched 24 3 x 64 
lattices, with a pion mass of m K = 227.6(6) MeV, a kaon mass of = 564(2) MeV, and pion 
momentum p % = ^fnnjL = HOy^MeV where n is the number of twists (0, 1, 2, or 3). It was 
necessary to reduce to the time separation between the kaon and the two pion sources to 20, 24, 
28, and 32 units. An error weighted average of the results from all four separations produced an 
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Table 2: Values of the two pion energies {E nn ) extracted from the effective mass plots. These are compared 
to the value one would expect in the absence of interactions between the two pions, 2^Jm\ + p 2 , as a sanity 
check. 





l^m^ + p 2 (lattice units) 


E nn (lattice units) 


2^/ml+p 2 (MeV) 


Enn (MeV) 


p = 


0.3501 


0.3509(13) 


455 


456(2) 


p = n/L 


0.4372 


0.4447(19) 


568 


578(2) 


p = V2n/L 


0.5096 


0.5316(50) 


662 


691(6) 


p = V3n/L 


0.5729 


0.577(20) 


745 


750(27) 



Table 3: Lattice matrix elements for the different time separations between the kaon and the two pions, 
At =20, 24, 28, and 32. The error weighted average of the results from these four separations is shown in 
the final column. Notice that its error is generally significantly smaller than the corresponding errors in the 
other columns. 



Matrix Elements with (27,1) Operator 




At = 20 


At = 24 


At = 28 


At = 32 


Error Weighted Ave. 


p = 


0.00529(29) 


0.00540(28) 


0.00490(26) 


0.00550(34) 


0.00526(17) 


p = n/L 


0.00882(55) 


0.00923(59) 


0.00813(57) 


0.00924(68) 


0.00884(37) 


p = V2n/L 


0.0140(11) 


0.0150(16) 


0.0137(17) 


0.0157(19) 


0.01449(99) 


p = V37T/L 


0.0249(47) 


0.0278(62) 


0.0191(49) 


0.0262(91) 


0.0240(38) 


Matrix Elements with (8,8) Operator 




At = 20 


At = 24 


At = 28 


At = 32 


Error Weighted Ave. 


p = 


0.0864(34) 


0.0917(37) 


0.0863(40) 


0.0947(45) 


0.0895(28) 


p = n/L 


0.0898(40) 


0.0918(41) 


0.0861(51) 


0.0951(57) 


0.0906(31) 


p = \/2n/L 


0.0895(69) 


0.0912(76) 


0.087(11) 


0.113(15) 


0.0931(61) 


p = yfin/L 


0.092(19) 


0.098(19) 


0.120(32) 


0.171(43) 


0.111(14) 


Matrix Elements with (8,8) mixed Operator 




At = 20 


At = 24 


At = 28 


At = 32 


Error Weighted Ave. 


p = 


0.297(11) 


0.315(12) 


0.295(13) 


0.324(15) 


0.3068(91) 


p = n/L 


0.319(14) 


0.327(14) 


0.302(17) 


0.333(19) 


0.320(11) 


p = V2n/L 


0.336(25) 


0.342(27) 


0.319(37) 


0.416(51) 


0.347(21) 


p = ^n/L 


0.377(68) 


0.385(74) 


0.45(14) 


0.66(15) 


0.437(51) 
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overall result with a smaller error. This error was of order 3-4% for pion momentum and Tt/L, 
7% for pion momentum y/2n/L, and 15% for pion momentum y/3n/L. 

Future plans include an expansion of this study (on the same lattices) to include a larger range 
of pion and kaon masses in order to perform an extrapolation to the physical pion and kaon mass, 
and to the physical (on-shell) pion momentum using the range of momenta obtained from the differ- 



ent twists. The Lellouch-Luscher factor [ |12| ] will be used to convert the lattice matrix elements to 
physical decay amplitudes. The calculation will also be performed on the RBC/UKQCD 32 3 x 64, 
L s = 32 lattices with 2+1 flavors of DWF that are currently being generated, with a pion mass of 
m n < 250 MeV, in order to obtain results in the full dynamical theory which can be compared with 
the results of the quenched approximation. 
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